Lifetime and Coherence of Two-Level Defects in a Josephson Junction 
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We measure the lifetime (Ti) and coherence (T2) of two-level defect states (TLSs) in the insulating 
barrier of a Josephson phase qubit and compare to the interaction strength between the two systems. 
We find for the average decay times a power law dependence on the corresponding interaction 
strengths, whereas for the average coherence times we find an optimum at intermediate coupling 
strengths. We explain both the lifetime and the coherence results using the standard TLS model, 
including dipole radiation by phonons and anti-correlated dependence of the energy parameters on 
environmental fluctuations. 



Two-level defects in amorphous insulators are of fun- 
damental interest due to their impact on many low tem- 
perature properties, such as the heat conductivity [1] and 
the generation of 1// noise [2, 3]. On the practical side, 
these effects limit the operation of solid state devices, 
for example amplifiers [4] and CCD detectors [5], and 
increase the dielectric loss of insulators [6] . 

Recently, a new type of solid state device has emerged 
in which quantum coherence is maintained over a large 
distance. In particular, superconducting Josephson 
qubits allow one to study quantum coherence at the 
macroscopic level [7]. Two-level defects in their amor- 
phous oxide tunnel barriers (usually made of A10 x (x ss 
1) [8]) have been found to limit the performance of these 
devices [6, 9, 10]. In the phase qubit, it was found that 
at certain biases the qubit is strongly coupled to spuri- 
ous two-level states (TLSs) which result in free oscilla- 
tions with the qubit and effectively reduce its coherence 
[11]. Similar effects have been observed in the flux qubit 
as well [12], although these are less common due to its 
smaller junction. Supercondcuting qubits hold promising 
features for the implementation of quantum information 
processing devices, however a significant improvement in 
defect density is required for future progress. 

These defects are thought to arise from charge fiuctu- 
ators in the insulating material of the junction, presum- 
ably O-H bonds [6] . Measurements on dielectrics at high 
temperature and power combined with measurements on 
the phase qubit at low temperature strongly support a 
two-level model for these fluctuators, emerging from tun- 
neling between two configurational states of the charge 
inside the junction. The notion of coupling of the phase 
qubit to a strongly anharmonic microscopic system was 
further fortified through careful analysis of the multilevel 
spectrum of the phase qubit near resonance with a defect 
[13]. Neeley et al. have demonstrated coherent control 
over a single defect, and characterized its coherence and 
relaxation time [14]. Other attempts are being made to 
reduce the impact of defects by improving junction ma- 
terials using epitaxial growth [15] or by fabricating di- 
electric free junctions [16]. 
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Figure 1: Frequency-domain and time-domain signatures of 
TLSs in qubit measurement, (a) The probability of the qubit 
being excited (Pi, normalized) after a long microwave tone 
at different frequencies and biases. On top of the expected 
smooth change of qubit frequency with bias, we observe ran- 
domly scattered splittings due to coupling to TLSs. (b) Pi 
after qubit excitation followed by a "free" evolution at different 
bias values. Pi oscillates near resonance with TLSs, consis- 
tent with the position of splittings in the spectrum. Upper 
and lower insets: control sequences used to produce (a) and 
(b) respectively. |0, g) and |l,g) stand for states where the 
qubit is at its ground and excited state respectively while the 
TLS is at its ground state. Dashed inset: level diagram for 
the combined qubit-TLS system on resonance. 



Several mechanisms have been proposed for relaxation 
and dephasing of the dielectric defects themselves. En- 



ergy relaxation is caused by coupling to phonon states, 
while dephasing could be caused by spectral diffusion 
[1]. However, to date, only limited measurements were 
carried out on TLSs to characterize these processes at 
the single defect level [17]. Such measurements could be 
used to better understand the nature of the defects and 
their decay mechanisms, and possibly engineer long-lived 
quantum memories in future devices. In [17], the coher- 
ence times of several TLSs were measured spectroscopi- 
cally and were found to distribute as P ~ I/T2, where 
Ti is the spectroscopic coherence time. Neeley's method 
[14] of probing the TLS adds the capability of measur- 
ing the coherence time more accurately and also measure 
their lifetime separately. 

In this letter, we present a measurement of the decay 
of energy and coherence for a large ensemble of TLSs in 
a small area junction using the phase qubit. We find 
that on average, the energy relaxation time (Ti), fol- 
lows a power law dependence on the coupling param- 
eter to the phase qubit. The exponent of this power 
law is in fair agreement with what is expected from 
phonon radiation by a dipole (proportional to the cou- 
pling strength) inside the junction. The average dephas- 
ing time (T«£ = (I/T2 — l/2Ti) ) is coupling dependent 
as well, peaking at intermediate couplings. We interpret 
this optimum coupling to be caused by anti-correlated 
fluctuations in the physical parameters which determine 
the TLS energy. 

For small area junctions (~ 1/im 2 ), the typical mea- 
surement bandwidth allows us to detect and measure 
about 10 TLSs in a particular cooldown. Instead of using 
many different samples to acquire sufficient statistics, we 
use the fact that heating resets the TLS characteristics. 
The device is thermally anchored to the mixing cham- 
ber of a dilution refrigerator during measurement. We 
find that the TLS distribution is reset upon raising the 
temperature above 20 K and cooling down to the base 
temperature (10 mK). Some memory of the TLS distri- 
bution remains if the temperature is increased to only 
1.5 K [18]. We utilize this feature to produce a new set of 
TLSs and generate an ensemble. The data was taken over 
82 different TLSs, obtained from 8 different cooldowns. 

The initial identification of TLSs and their coupling pa- 
rameters are carried out as follows. First, the qubit spec- 
trum is swept over the bias [22] to locate the frequencies 
of TLSs from the positions of the avoided level crossing 
structures (see Fig. la). A complementary picture of 
the interacting qubit-TLS system in the time domain is 
shown in Fig. lb, where we excite the qubit with a short 
resonant pulse (7r-pulse) far away from any observable 
TLS and then apply a bias pulse of varying amplitude 
and length. As seen in the figure, for bias values where 
the qubit is resonant with a TLS we observe oscillations 
that have the same frequency as the splitting size in the 
spectrum. 

Following Neeley et al. [14], the characteristic energy 
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Figure 2: Representative Ti and T2 measurements of a TLS. 
(a) Ti measurement of the TLS, along with its experimen- 
tal sequence (inset). The qubit is first excited with a 7r-pulse, 
then brought into resonance with a TLS for a "swap time" (the 
time to fully transfer an excitation between the qubit and the 
TLS. It is found for each TLS by locating the first minimum 
in the oscillations in Fig. lb). After a free evolution of the 
TLS of time r, a swap gate is again applied, afterwhich the 
qubit excitation probability Pi is measured, (b) T2 measure- 
ment of a TLS, along with its experimental sequence (inset). 
This sequence is similar to a T\ measurement, only that su- 
perposition states are produced in the TLS and their decay 
is measured vs. time (Ramsey sequence [14]). The amplitude 
of the oscillations is proportional to the degree of coherence 
in the TLS. 



relaxation and decoherence time scales were extracted 
from T\ and Ramsey experiments on the TLS, with se- 
quences schematically represented in the insets of Fig. 
2. Figure 2a and 2b show representative T\ and T2 de- 
cay curves of the same TLS with characteristic times of 
187 ns and 148 ns respectively, obtained from a fit to a 
decaying exponent and an oscillatory decaying exponent. 
The size distribution of the observed splittings (see Fig. 
3a) follows the theoretical curve predicted by the stan- 
dard model for two- level defects and agrees with previous 
results on similar junctions (generated by measuring dif- 
ferent samples) [6] . The maximal splitting size is found to 
be 45 MHz. Theory predicts [6] that this maximal split- 
ting Smax depends on junction parameters and defect size 
according to S m&x = (2d/x) \Je 2 E w /2C , where x is the 
barrier thickness, d is the spatial size of the dipole, C 
is the junction capacitance and Eiq is the qubit energy. 
From the measured S max and known junction parame- 
ters we compute a dipole size d ~ 1 A. The minimal 
observable splitting size is ~ 3 MHz, and is mainly lim- 
ited by the coherence time of the qubit. In addition, we 
find the distribution of TLS energies (-E ge , the energy be- 
tween the ground state and excited state) to be constant 
throughout our qubit measurement bandwidth (see Fig. 
3b), consistent with theory. 

Although most of the T\ decay curves of the TLSs are 
similar in their shape (i.e. a simple exponential decay), 
their decay times range almost 3 orders of magnitude, 



from 12 ns to more than 6000 ns. Coherence times on the 
other hand range from 30 ns to only 150 ns (excluding a 
single anomalous TLS which will be discussed later). For 
comparison, when the qubit is biased far from any observ- 
able splitting, its lifetime is 270 ns, and its coherence time 
is 90 ns. TLS energy relaxation times at a given splitting 
are not random. We find that they are shorter for larger 
splittings (stronger interaction with the qubit) , although 
short lifetimes are measured for the smallest splittings as 
well. This trend is apparent in Fig. 3c where we plot 
average T\ values as a function of splitting. In this plot 
we divide the ensemble into groups of TLSs having sim- 
ilar splitting values, in a 7 MHz window size. The error 
bar represents the statistical spread of the data within 
this window [18]. We find the average values (Ti(5)), 
excluding two points [18], to be best fitted by a power 
law Ti a 5 Q , where a = -1.44 ± 0.15 [18]. Figure 3d 
(black circles) shows the processed T% data, obtained sim- 
ilarly from only 43 different TLSs [18]. In this case we 
observe a weak dependence on the coupling with a peak 
at S ~ 25 MHz. This feature is more pronounced in the 
dephasing time X^, represented by red triangles in Fig. 
3d. 

The T\ results can be understood within the standard 
TLS model. The excited state of the TLS involves a local 
deformation of the insulator. This deformation couples 
to phonon modes, leading to the decay of the TLS ex- 
citation. The expected lifetime for such a process [1], is 
given by 
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Figure 3: TLS survey results. (a) Frequency and (b) 
splitting-size distribution of 82 TLSs. The curve in (a) is 
the best-fitted log-normal distribution of the standard TLS 
model [6]. (c) Average Ti values (black circles) as a function 
of average splittings taken for TLSs inside a 7 MHz splitting 
window, and the best fitted power law in dashed gray. Av- 
erage Ti from stochastic simulation (arbitrary amplitude is 
shown) in blue diamonds. We attribute the deviation at the 
largest splittings to low statistics within these windows [18] 
(d) Average Ta (black circles) and T/> (red triangles) value as 
a function of average splitting size within a 7 MHz window. 
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where 7 is the deformation potential, v\ and vt are the 
speeds of sound for the longitudinal and transverse modes 
respectively, Ao is the energy splitting due to tunneling 
and p is the mass density. This is consistent with a power 
law dependence on S, since the interaction strength with 
the qubit satisfies S ex S ma ^A /E ge [6]. 

The interaction of a TLS with the qubit is that of an 
electric dipole with an electric field, and therefore de- 
pends on the dipole orientation [6] . This feature explains 
the large spread in the data at a given splitting: both 
large dipoles (large A /E ge ) perpendicular to the junc- 
tion's electric field and small dipoles aligned with the 
field can have the same S but different lifetimes. To more 
rigorously compare experiment with this theory, we sim- 
ulate an ensemble of TLSs with uniform distribution of 
dipole orientation and log distribution of dipole moment 
sizes [6] , from which we calculate the average lifetimes as 
a function of splitting size [18]. The simulation data (see 
Fig. 3c) yields an average exponent a = —1.63 [18], in 
agreement with our measurement. 

The magnitude of the times that we extract from the 
experiment at a given splitting can be compared to the 
expected values for defects inside an A10 x dielectric using 



Eq. 1. We approximate the deformation potential by 
7 = ^pv 2 AV [19], where v is the average speed of sound 
of the transverse and longitudinal modes and AV is the 
local difference in volumes. For AI2O3 values with the 



difference in volumes taken as AV 



where a ~ 1 A 



is the extracted dipole size, we get 7 saleV, consistent 
with defects in other dielectrics [3]. Since the dielectric 
layer thickness is much smaller than the relevant phonon 
wavelength, the speed of sound in Eq. 1 is set by the 
aluminum layers of the electrodes. Using the speed of 
sound for thin aluminum films [20], we get Ti(S max ) ~ 
30 ns which is very similar to what we measure. A more 
specific estimation should take into account the size of 
the junction and the layered structure. 

Assuming the dephasing process is caused by fluctu- 
ations in energy, we note that the maximum observed 
in Fig. 3d can be explained by an anti-correlated de- 
pendence of the charging energy and tunneling energy 
on fluctuations in the TLS environment. According to 
the TLS model. E a 



v/A 



1 - 1 - Aq where A is the en- 



ergy difference between the bare states of two spatial 
configurations \L) and \R) and A is the tunneling in- 
teraction energy. Both A and Ao are dependent on a set 
of environmental parameters P , which fluctuate in time. 
As is standard for the TLS model, we assume a linear 
sensitivity for A on P and an exponential sensitivity 



for A : A (7*) - N ie -^ p *^, A(P) = N 2 ^Pi/P u , 
with overall dimensional normalization constants Ni and 
N 2 and parameter specific constants Pq and Pi. The 
resulting fluctuations in energy, to first order in fluc- 
tuations SP in these parameters, are given by 6E Ee = 
-gr— ^2 SPi (A 2 /Pii — Aq/Poj)- This expansion becomes 
interesting for the situation where Ao ~ A as there is a 
possibility for the contribution in 5E gc to pairwise can- 
cel. Since A = S {E gc /S max cos 77), where 77 is the dipole 
orientation relative to the junction's electric field, we ex- 
pect to find such a cancellation at a particular splitting 
S. Note that the dependence on cos 77 smears this some- 
what but we still expect a significant effect, as is observed 
in Fig. 3d. 

As seen in the figures, the power-law describing Ti(S) 
cannot explain all the measured TLSs. We find that three 
out of 82 TLSs with large splittings (37 MHz, 41 MHz 
and 45 MHz) have much longer lifetimes than expected 
(220ns, 243ns and 476ns - respectively). In addition, 
one TLS out of 41 has much longer coherence than all the 
others (about a factor of 6 longer than the longest T 2 of 
all the others), associated with a splitting of size 30 MHz. 
Other anomalies we discovered are related to the stability 
of a particular TLS in time. We find that the energy 
.Ege of some TLSs (about 5 %) changes spontaneously at [8 



varying time scales, from seconds to days. All the rest 
were remarkably stable [18]. 

Some of these changing TLSs have long lifetime (a few 
microseconds), which is consistent with the power-law 
trend we discussed above. Furthermore, we also measure 
a few representative TLSs as a function of temperature. 
We find no significant change in T\ and T 2 below 100 mK, 
consistent with the expected [1] tanh(_E g0 /2fcBT) depen- 
dence. We also find that the instability of some TLSs 
increases at elevated temperatues (i.e, the change in TLS 
energy becomes more frequent). We conclude that some 
of the TLSs we measure have a different nature, perhaps 
related to their internal structure or position inside the 
junction. 

In conclusion, the energy decay and dephasing times 
of two-level defects in an A10 x barrier of a Josephson 
junction are measured as a function of the coupling pa- 
rameter with the phase qubit. The lifetimes vary sub- 
stantially in our range of splittings, and agree with the 
theoretically predicted phonon radiative loss, which is 
dipole size dependent. The dephasing times show an ex- 



trcmum at intermediate couplings, which we attribute to 
an anti-correlated dependence on fluctuations in the en- 
vironmental parameters which set the TLS energy. Such 
a dependence may distinguish between different theoret- 
ical models for TLSs. Our results demonstrate the power 
of the phase qubit as a dynamical coupling element to mi- 
croscopic systems at the single microwave photon level. 
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EXPERIMENTAL DATA 

The full data measured for 82 TLSs is presented in Ta- 
ble I and the Ti(S), T 2 {S) and T (5) values are plotted 
in figure 1. Figure la contains 82 data points while lb 
contains only 42. For about 50 % of the cases the de- 
phasing time could not be determined, due to both low 
visibility and short dephasing time. This happens mostly 
for small splittings, which do not lie in the range of points 
presented in Fig. 3d of the main paper. For points that 
lie within this range and are omitted from the analysis 
we separately checked that the shorter coherence time 
does not affect the trend. The number of points in Fig. 
lc is further reduced because for some T\ limited TLSs 
we measured T 2 which is slightly longer than 2T\ due to 
measurement error. T^ is excluded from the figure and 
table for these cases (three TLSs). 
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Figure 1: (a) Measured T\ vs. splitting, (b) T2 vs. splitting 
and (c) T<p vs. splitting for all the measured TLSs. The spread 
in the T\ data at a particular splitting value results from 
the random distribution of TLS orientation in the junction. 
It is apparent that the maximal lifetime shortens at larger 
splittings, consistent with dipole radiation. The spread of 
the T2 data appears independent of the splitting, however the 
average values show some dependence which becomes more 
prominent in the dephasing times. 

The errors in Fig. 3b and and Fig. 3c in the main 



paper represent the statistical spread of the data within 
a 7 MHz window. They are calculated by normalizing 
the standard deviation by vN, where N is the number 
of points within the window. 
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Figure 2: Qubit spectrum as a function of bias, before and 
after warmup to 1.5 K. Upper: two sections of the spectrum 
before warmup. Lower: the same sections, taken after heating 
the sample to 1.5 K and cooling down back to lOmK. Circles 
with the same color indicate a splitting that we attribute to 
the same TLS. Both the position in frequency and splitting 
size of the TLSs are similar, indicating that the TLSs are not 
fully reset at this temperature. 



We would like to point out that two points in Fig. 3c 
in the main paper have been excluded from the fit to a 
power law. These are the points with the largest aver- 
age splittings (38 MHz and 46 MHz), where the statistics 
within each window is low (3 data points and 2 data 
points respectively). If we include these points in the 
fit we obtain an exponent a — —0.93. We believe that 
the data points at the largest splittings may result from 
anomalous TLSs. 

As pointed out in the main paper, some of the TLSs are 
changing in time. This is characterized by a change in the 
TLSs' transition energy, causing them to disappear from 
time to time (that is, to step out of our measurement 
bandwidth) or appear at a slightly different energy. An 
example of this phenomenon can be seen in Fig. 1 of the 
main paper. The leftmost TLS appears in the spectrum 
at larger bias value (that is smaller transition energy) 
than in the time-domain sweep. These two measurements 
were taken at an interval of one day. 

New sets of TLSs are produced from the same device by 



warming up to 20 K. We believe that some TLSs are not 
fully reset after warming up to only 1.5 K. Some splittings 
in the spectrum are similar in their frequency and split- 
ting size to those before a partial warmup, as indicated 
in Fig. 2. Both the reset of TLSs at high temperatures, 
and the fact that some of the TLSs are changing in time 
indicate the true nature of the TLS as an approximation 
of a multilevel state, resulting from a multi-well energy 
structure. 



STOCHASTIC SIMULATION 



(A oc z, where z is the effective distance between posi- 
tion states of the TLS inside the junction) and a tunnel- 
ing energy Ao from a log distribution (In Ao oc z) . These 
distributions are consistent with the TLS model, which 
assumes linear sensitivity of A on z and exponential sen- 
sitivity of Ao on z. We excluded TLSs having smaller 
couplings than we can measure. For each TLS we com- 
pute the lifetime according to Eq. 1 in the main paper 
(Ti(sin(0)) = a/sin(6>) 2 , where 9 = arctan(A /A) and 
a is some constant). The simulation data points were 
then averaged over a 7 MHz window size, as done for the 
experimental data. 



We reproduced the lifetime distribution of an ensem- 
ble of TLSs according to the TLS model in the following 
way. For each TLS we assume a uniform distribution of 
dipole orientation (cos rj distributes uniformly, where r\ 
is angle relative to the electric field inside the junction), 
a potential asymmetry A from a uniform distribution 



The resulting (T\(S)) behavior (see Fig. 3c in the main 
paper, blue diamonds) resembles a sum of two power 
laws. At smaller splittings S < S max /2, the points fit a 
power law with an exponent ol\ sa —1.9, while for larger 
splittings they fit a power law with an exponent a 2 ~ — 1- 
A similar trend is observed in our data as well. 
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5 


11.62 


16.1 


70 


100 


350 


5 


10.95 


16.1 


600 


- 


- 


1 



/ge [GHz] 


S [MHz] 


Ti [ns] 


T 2 [ns] 


T* [ns] 


CD No. 


12.642 


16.4 


1080 


70 


72 


2 


11.85 


16.4 


37 


- 


- 


7 


12.41 


17.5 


165 


110 


165 


5 


12.37 


17.5 


150 


- 


- 


7 


12.542 


18.0 


106 


81 


131 


8 


11.42 


18.5 


91 


144 


690 


4 


11.915 


18.9 


208 


- 


- 


2 


12.47 


19.2 


72 


36 


48 


6 


11.21 


20.0 


127 


100 


165 


5 


12.73 


20.8 


41 


- 


- 


7 


11.86 


21.7 


156 


120 


195 


2 


12.15 


22.2 


85 


- 


- 


7 


11.72 


23.8 


60 


120 


- 


2 


12.66 


23.8 


32 


- 


- 


7 


12.44 


25.0 


44 


- 


- 


6 


11.86 


25.0 


355 


- 


- 


1 


11.61 


25.6 


53 


100 


1767 


3 


12.34 


25.6 


32 


- 


- 


6 


11.73 


25.6 


224 


48 


54 


7 


13.04 


26.3 


217 


50 


57 


6 


12.675 


27.0 


110 


- 


- 


1 


11.77 


27.8 


62 


130 


2687 


4 


12.2 


28.0 


50 


- 


- 


1 


12.613 


29.4 


56 


52 


97 


2 


11.22 


29.4 


89 


41 


53 


2 


12.67 


29.4 


193 


50 


57 


3 


12.147 


30.0 


150 


212 


723 


8 


12.34 


30.3 


82 


58 


90 


4 


11.89 


31.3 


31 


- 


- 


6 


10.91 


31.3 


300 


120 


150 


7 


12.78 


34.5 


125 


114 


210 


5 


12.12 


34.5 


48 


52 


113 


6 


11.59 


35.7 


26 


27 


56 


2 


11.17 


37.0 


243 


68 


79 


6 


11.772 


41.0 


220 


156 


242 


1 


12.05 


45.5 


476 


950 


- 


4 


11.99 


47.0 


60 


- 


- 


1 



Table I: Full measurement data of TLSs: TLS energy (/ ge ), splitting (S), lifetime (Ti), coherence time (T 2 ) calculated dephasing 
time (T<t>) and cooldown number. The last column indicates which TLSs were measured on the same cooldown (i.e., temperature 
was not raised to more than 20mK between measurements of TLSs belonging to the same cooldown). 



